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By means of two-dimensional contact dynamics simulations, we analyze the vibrational dynamics 
of a confined granular layer in response to harmonic forcing. We use irregular polygonal grains 
allowing for strong variability of solid fraction. The system involves a jammed state separating 
passive (loading) and active (unloading) states. We show that an approximate expression of the 
packing resistance force as a function of the displacement of the free retaining wall from the jamming 
position provides a good description of the dynamics. We study in detail the scaling of displacements 
and velocities with loading parameters. In particular, we find that, for a wide range of frequencies, 
the data collapse by scaling the displacements with the inverse square of frequency, the inverse of the 
force amplitude and the square of gravity. Interestingly, compaction occurs during the extension of 
the packing, followed by decompaction in the contraction phase. We show that the mean compaction 
rate increases linearly with frequency up to a characteristic frequency and then it declines in inverse 
proportion to frequency. The characteristic frequency is interpreted in terms of the time required 
for the relaxation of the packing through collective grain rearrangements between two equilibrium 
states. 

PACS numbers: 



I. INTRODUCTION 

Depending on the frequency and amplitude of accel- 
erations, vibrated granular materials give rise to vari- 
ous phenomena such as compaction, ([II 0) convective 
flow(|3, 01 I3)j size segregation and standing wave pat- 
terns at the free surface (J., ^ JIj). Particle rearrange- 
ments induced by vibrations lead to lower shear strength 
and larger flowability. In the full fluidization regime, 
there are no permanent contacts between particles and 
the system behaves as a dissipative gas Q- Particle bed 
reactors are sometimes fluidized by this method instead 
of upward gas flow (Q). When particle accelerations 
remain below the gravitational acceleration, the system 
keeps its static nature and the vibrational energy propa- 
gates through a rather compact network of interparticle 
contacts. This leads to enhanced bulk flow in hoppers 
and chutes ilfl, ilTj . 

On the other hand, vibrations at high frequency and 
low amplitude lead to slow (logarithmic) decay of the 
pore space as a function of time [5J- Efficient vibro- 
compaction of dry and wet granular materials is a cru- 
cial issue in numerous applications such as the casting of 
fresh concrete. The tamping operation on railway bal- 
last is another example where the vibrations of tamping 
bars are used to restore the initial geometry of the track 
distorted as a result of ballast settlement (|12. [l3. [l4l|). 
The maintenance cost becomes crucial with the increase 
of commercial speed. 

We may distinguish two methods for inducing vibra- 
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tional dynamics: 1) by imposed cyclic displacements of a 
wall or the container (shaking); 2) by cyclic modulation 
of a confining stress. The first method has been used in 
most experiments on granular beds ( 15, 16, 17, 18, 13, 
I20I'). In this case, the control parameters are the ampli- 
tude a and the frequency v of the vibrations correspond- 
ing to a maximal acceleration aw^ where uj = 2itv. When 
a material is moulded inside a closed box, the vibrations 
should rather be induced by varying a confining force, 
e.g. a force acting on a wall. Then, the amplitude of 
displacements is a function of the forcing frequency, and 
the level of particle accelerations depends on both the 
applied cyclic force and the reaction force of the pack- 
ing. In any case, an efficient compaction process requires 
periods of release of the packing so that the grains can 
move with respect to their neighbors. 

In this paper, we explore such a system where a har- 
monic force / is exerted on a lateral wall of a closed box, 
all other walls remaining immobile. The force / is var- 
ied between zero and a maximum value fmax- During a 
period, / is large enough to equilibrate the packing reac- 
tion force except for a short laps of time when / declines 
to zero. Then, the packing can flow under the action 
of its own weight, pushing the retaining wall away. We 
are interested here in the evolution of the packing in the 
course of harmonic loading and its scaling with loading 
parameters (frequency, force maximum). 

We used numerical simulations by the contact dynam- 
ics approach as a discrete element method (DEM) in a 
two-dimensional geometry with a small number of par- 
ticles (I21II23). Each simulation is repeated for several 
independent configurations and the results are analyzed 
in terms of ensemble average behaviors. The particles 
are rigid and polygon-shaped. We focus on the displace- 



merits of the free retaining wall and the compaction of 
the packing. Most results presented below concern the 
short-time behavior where the solid fraction increases lin- 
early with time. The frequency is varied from 1 to 60 Hz 
and its influence is analyzed by considering characteristic 
times involved in the loading and unloading intervals of 
time. We first introduce the numerical procedures. Then, 
we present the main findings concerning the passive and 
active dynamics, the evolution of the solid fraction and 
scaling with the loading parameters. 



II. NUMERICAL PROCEDURES 

The simulations were carried out by means of the con- 
tact dynamics (CD) method with irregular polygonal 
particles (' |2ll ^^). The CD method is based on im- 
plicit time integration of the equations of motion and a 
nonsmooth formulation of mutual exclusion and dry fric- 
tion between particles. This method requires no elastic 
repulsive potential and no smoothing of the Coulomb fric- 
tion law for the determination of forces. For this reason, 
the simulations can be performed with large time steps 
compared to molecular dynamics simulations. We used 
LMGC90 which is a multipurpose software developed in 
our laboratory, capable of modeling a collection of de- 
formable or undeformable particles of various shapes by 
different algorithms (|23|)- 

The samples are composed of irregular pentagons, 
hexagons, and heptagons of three different diameters: 
50% of diameter d„iin = 2.5 cm, 34% of diameter 3.75 cm, 
16% of diameter dmax = 5 cm; see fig. ^ The particles 
are initially placed on a square network in a rectangular 
box and compressed by downward motion of the upper 
wall (wall C in fig. QJ at zero gravity. Then, the gravity 
is set to g and the upper wall is raised 1 cm and fixed. 
The right wall (wall D in fig. ^ is allowed to move hor- 
izontally (x direction) and subjected to a driving force: 
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where fmax and /mm sue the largest and lowest compres- 
sive (positive) forces acting on the wall. 

If fmin is above the (gravitational) force exerted by the 
grains on the free wall, / will be large enough to prevent 
the wall from backward motion during the whole cycle. 
In other words, the granular material is in "passive state" 
in the sense of Rankine's states and the major principal 
stress direction is horizontal (f24j). In this limit, no ex- 
tension will occur following the initial contraction. On 
the other hand, if fmax is below the force exerted by the 
grains, / will never be large enough to prevent the ex- 
tension of the packing. This corresponds to the "active 
state" where the major principal stress direction remains 
vertical. In all other cases, both contraction and exten- 
sion occur during each period, and the displacement Ax 
of the free wall will be controlled by fmin- In the sim- 
ulations reported below, we set fmin = 0. This ensures 




FIG. 1: The geometry of the packing. 



the largest possible displacement of the wall in the active 
state. We used four different values of fmax ranging from 
5 .10^ N to 2 .10"^. 

The simulations were carried out with Np — 95 grains 
in the box and each simulation was repeated with seven 
independent grain configurations. The mean behavior for 
each set of parameters is obtained by ensemble averaging 
over seven independent data sets. Larger samples can be 
simulated, but that requires much more computational 
effort for a parametric study over many cycles. Thus, 
our system represents rather a thin granular layer. The 
coefficient of friction between the grains and with the 
horizontal walls was fixed to 0.4, but it was at the ver- 
tical walls. With a time step equal to 2.5 10~* s we could 
perform high-quality simulations in which the largest cu- 
mulative error on grain positions was bellow 1%. 



III. ACTIVE AND PASSIVE DYNAMICS 

We first consider the motion x{t) of the free wall (wall 
D in Fig. ^ which reflects the dynamics of the grains in 
the cell in response to harmonic forcing. Figure |21 shows 
x{t) (averaged over seven independent simulations) for 
frequency v — 5 Hz over a time interval At = 1 s. We 
distinguish a fast initial contraction (t < 0.1 s) followed 
by slow contraction (decreasing x) over four periods. The 
initial contraction is a consequence of the gap left be- 
tween the free surface of the packing and the upper wall. 
This initial volume change is almost independent of fre- 
quency. The subsequent periodic motion of the wall takes 
place around this conflned state and will be at the focus 
of this paper. 

A zoom on a single period is shown in Fig. |31 The pe- 
riod begins at the jamming position x — Xj correspond- 
ing to the jamming position reached in the preceding pe- 
riod. The motion of the wall begins (point a in Fig. O 
only when the applied force / declines near to its min- 
imum fmin = 0. The maximum displacement Axmax 
occurs at a later time dtr (point b). From a to b, the 
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FIG. 2: The evolution of the displacement x of the free wall 
(up) in response to harmonic loading (down). 
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FIG. 3: A zoom on a single period; see Fig. |21 



force exerted by the packing on the free wall is above the 
applied force, so that the wall moves backward (exten- 
sion). In this phase, the packing is in an active state. The 
inverse situation prevails from b to c where the grains are 
pushed towards the box (contraction). Then, the pack- 
ing is in a passive state. The new jamming position x'- 
is below the jamming position Xj reached at the end of 
the preceding period. The difference Xj — x'- represents 
the net compaction of the packing over one period. For 
a given frequency v, the phase difference Str is the same 
for all periods. The displacement amplitude Axmax is a 
function of fmax and v, as we shall see below. 

The motion of the free wall is governed by the equation 
of dynamics. 



f - .fg = mx, 



(2) 



where fg is the horizontal force exerted by the packing 
on the wall (fig^. Figure 0] displays fg as the function of 
time for fmax — 2 10'* N. We see that fg follows closely 
the variations of /. In particular, in the jammed state 
we have f = fg so that i = in this state. This means 
that, in its most general form, fg is a function of /. 
Figure EJa) shows fg vs. x over four periods. In the 



\ 1 




r 




-fg 


iJ 




--f 


M 


\ 




- 


\ 'f 

- \ ' f — 

\ 'J 


1 ^Hi*cr^ 1 



<u 
o 

^3 



07400 0.425 0.450 0.475 0.500 

t(s) 

FIG. 4: The force fg exerted by the grains and the driving 
force / on the free wall as a function of time i. 
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FIG. 5: Force fg exerted by the grains on the free wall as a 
function of displacement x (a) and the velocity v (b). ( 



active phase, fg grows slightly with x. In the passive 
phase, it grows faster and almost linearly as x decreases. 
The vertical line corresponds to the jammed state where 
fg decreases with f at x = Xj. We also clearly observe in 
Fig. EJa) two transients : 1) unjamming and the onset 
of the active state, 2) jamming from the passive state. 





FIG. 6: Particle displacements over one period. 



It is remarkable that, although Xj decreases at the end 
of each period, the dynamics remains self-similar up to a 
translation along displacement coordinates. 

Figure E^b) displays fg as a function of the velocity 
V = X. We again observe the passive (v < 0) and active 
(v > 0) states together with the jamming and unjamming 
transients before and after the jammed state {v = and 
X = Xj). The data from all periods follow the same 
variations except for the jamming transient where a slight 
decrease of the maximum negative velocity Vmax can be 
noticed in each period. 

Although we focus here on the average dynamics of the 
packing, i.e. the displacements of the free wall, it is im- 
portant to note that the grain velocity field is not a simple 
oscillation around an average position. The grains un- 
dergo a clockwise convective motion in the cell as shown 
in Fig. El On the other hand, the contact forces evolve 
between a fully jammed state, where nearly horizontal 
force chains dominate (Fig. EJa)), and the active state, 
where nearly vertical gravity-induced chains can be ob- 
served (Fig. Ulh)). 



IV. A PHENOMENOLOGICAL MODEL 

To predict the motion of the free wall from equation 
12 we need to express the force fg as a function of x and 
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FIG. 7: Normal forces in the passive (a) and active (b) states. 
Line thickness is proportional to the force. 



peaks in each period correspond to the maxima of the 
largest velocities in absolute value shown in Fig. EJb). 
We see that / < 0.02, implying that fg should not depend 
crucially on v. Let us note that the plot of fg vs v in 
FigEIb) does not represent the explicit dependence of fg 
on w; it is a consequence of the equation of dynamics and, 
as we shall see below, it can be reproduced by assuming 
no dependence of fg on u. 

We now introduce a simple phenomenological model 



It is obvious that in the jammed state at x = Xj, the in which the expression of fg as a function of x is ex- 
tracted from the numerical data plotted in Fig. |3fa). As 
shown in Fig. |51 two distinct fitting forms are to be con- 
sidered for the active and passive states. Ignoring the 
jamming and unjamming short transients, an exponen- 
tial form provides a nice fit for the active branch whereas 
a linear fit seems fairly good in the passive state. Hence, 
with a good approximation we can write 



force fg is a reaction force balancing exactly the driving 
force / so that v — 0. On the other hand, the inertia 
effects are small compared to static forces. To show this, 
we may use a dimensionless number I defined by (,2^J): 



I ^i 



(3) 



where i = x/x is the deformation rate, m is the total 
mass, and p is the average pressure. The evolution of 
/ is plotted as a function of time in Fig. |S1 The two 



/.= 



a + (3 e'' (^-^^) 
a' + (3' {1 + k' {x 



active, 



Xj)} passive 



(4) 



0.020 



0.015 



i-H 0.010 



0.005 - 



OooqI i — LJ_U U_±A 1 J I V — ULLA 1 

0.0 0.2 0.4 0.6 0.8 1.0 
t(s) 

FIG. 8; The evolution of the inertia number / over five peri- 
ods. 
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FIG. 9: Variation of the packing reaction force fg with dis- 
placement X over one period (full line) and an approximate 
fitting form (dashed line). 
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The constant forces /i , /2 and /a correspond to the values 
of fg at the unjamming transient, the jamming transient 
and the point of transition from active to passive states, 
respectively (see Fig. ^. Clearly, because of the action 
of gravity and jamming transition, we have /i > fmin 

and /2 < fmax- 

We substitute this expression 0] in equation [3 and we 
solve for x. Analytical solution can be obtained for the 
passive linear part. An approximate solution can be given 
also for the active part by expending the exponential 
function to leading order. Figure E| shows the evolu- 
tion of the position x for one period together with the 
solution of the model. 

The parameters fc, k' are adjusted in order to get the 
best fit for the plot. The continuity of the fit at transition 
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FIG. 10: Displacement x of the free wall as a function of time 
(full line) and analytical fit from the phenomenological model 
(dashed line). 
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FIG. 11: Force fg exerted by the grains on the free wall versus 
velocity v over one period (full line) and analytical fit from 
the phenomenological model (dashed line). 
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FIG. 12: Phase space trajectory over one period (full line) 
and analytical fit from the phenomenological model (dashed 
line) . 
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FIG. 13: Evolution of the solid fraction Ap from the initial 
state as a function of time over several periods. 
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FIG. 14: Evolution of the jamming position Xj from the posi- 
tion x'^' reached at t = t'^' = 0.1 s for four different frequen- 
cies. 



between passive and active states is ensured by the very 
choice of the coefficients according to Eq. Although 
we did not take into account the transients, the analytical 
plot fits correctly the data. 

Figure im displavs fg vs v for one period, together with 
the analytical fit obtained as solution to Eq. [21 given the 
expression 0] of fg as a function of x. Again, excluding 
jamming and unjamming transients, the analytical solu- 
tion provides a fairly good approximation for the simu- 
lation data. Fig. 1121 shows the trajectory of the motion 
in the phase space (i, x) for one period, both from di- 
rect data and the model. The fit is globally acceptable 
although the velocity is under-estimated in the passive 
state. 

The model parameters k and k' remain nearly the same 
over all periods. This means that the dynamics at short 
times {At < 1 s) is weakly dependent on the solid frac- 
tion. The parameters k and k' change, however, with 
loading parameters (i/, fmax, etc) unless the displace- 
ments and the forces /i, /2 and /a are scaled with these 
parameters. This point will be discussed in detail below. 



V. COMPACTION 

In order to evaluate the solid fraction p, we consider a 
control volume enclosing a portion of the packing inside 
the simulation cell. This volume does not include the 
initial gap between the top of the packing and the upper 
wall. The initial value of the solid fraction is 0.75 and, 
since the grains are angular-shaped, its variations Ap 
from the initial state are large. 

Figurein^shows the evolution of Ap for several periods. 
We observe an initial compaction of 3% occurring in 0.1 
s. The subsequent evolution of the solid fraction takes 
place in a more compact state with a small increase in 
each period. 

We use po = 0.77 at the end of the first period as 
the reference value for solid fraction. The relative com- 



paction of the packing is given by Ap/po. The com- 
paction rate rj over several periods and for a total time 
interval At is 



1 Ap 



(6) 



Fig. 1 141 shows the jamming position Xj as a function of 
time for different frequencies for At < Is. At such short 
times, it can be assumed, with a good approximation, 
that the solid fraction declines linearly in time. Gener- 
ally, the behavior slows down logarithmically at longer 
times 16l|. This means that at short times, in which we 
are interested in this paper, the compaction rate is nearly 
constant, and we have 



Api 
Po 



(7) 



5 Hz 



where Api is the compaction per period. For u 
and fjnax = 2 10"* N, we have ry ~ 0.009 8"^ 

Interestingly, compaction occurs in the active state, 
i.e. during the extension of the packing, and not during 
contraction! This is shown in Fig. 1151 where the variation 
Ap of the solid fraction is plotted as a function of x. The 
solid fraction increases during extension (increasing x) 
and decreases during contraction (decreasing x). 

Compaction upon the reversal of the direction of shear- 
ing is a well-known property of granular media |2g . Low- 
amplitude cyclic shearing leads to cumulative compaction 
of a granular material. At larger amplitudes, the com- 
paction is followed by decompaction (dilation) and no net 
compaction can be observed over a full cycle. The sit- 
uation is slightly different in our system in the presence 
of a jammed state. Compaction is a consequence of un- 
jamming and it is pursued during the whole active state. 
Decompaction takes place in the passive state, but it is 
cut short by fast jamming. The outcome of a full cycle 
is thus a net compaction of the packing. 
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FIG. 15: Variation Ap of the solid fraction from the initial 
state as a function of the displacement x of the free wall. 



VI. SCALING WITH LOADING PARAMETERS 

In the last three sections, we analyzed the vibrational 
dynamics and compaction for a single frequency u = 5 
Hz. Similar simulations were performed for several fre- 
quencies ranging from 1 Hz to 60 Hz. Up to a change in 
time and length scales, all simulations yield similar re- 
sults both for dynamics and compaction independently 
of the applied frequency. This can be seen, for example, 
in Fig. llBf a) where the phase space trajectory is shown 
for t/ = 5 Hz and ly — 10 Hz. Fig. llGr b) shows that the 
data from both simulations collapse nicely on the same 
curve by simply scaling the displacements Ax by iy~^ and 
the velocities v by v^^ . 

This scaling is suggested by a dimensional analysis of 
the average dynamics of the packing. The frequency ly 
sets the time scale r = ly^^. Force scales are set by the 
largest driving force fmax in the passive state and the 
grain weights mg as well as the smallest driving force 
fmin in the active state. Hence, dimensionally, for fixed 
values of mg, fmin and fmax, all displacements are ex- 
pected to scale with j/~^ and all velocities with i^^^. To 
check directly this scaling, in Fig. El we have plotted 
the maximum displacement Axmax in the active state 
and the maximum velocity Vmax in the passive state as 
a function of v. The corresponding fits by z/~^ and v~^ 
are excellent. 

The influence of loading force parameters mg, fmin 
and fmax should be analyzed separately for each regime. 
In the passive state, fmax is the dominant force and it 
is exactly balanced by fg in jamming transition. On the 
other hand, in the active state, mg is the dominant force 
as / remains small compared to mg in this state. The 
maximum displacement Axmax at transition from active 
to passive state is determined in a subtle way by both 
fmax and mg. If gravity were the only driving force in 
the active state, Axmax would simply scale with gv'^ 
independently of fmax- However, our data show that 
Axmax varies as fmax, Fig. El A plausible dimensional 
interpretation is to assume that Axmax is controlled by 
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FIG. 16: Phase space trajectories for two frequencies with- 
out scaling (a) and with scaling (b) of the displacements and 
velocities. 



the ratio trig/ fmax representing the relative importance 
of the gravitational to loading forces. Then, the following 
simple expression can be proposed for the scaling with 
loading forces: 



Ax^ 



C 



mg 

/ma a 



m^ 



(8) 



where C is a dimensionless pref actor. This equation in- 
cludes the correct scaling of Axmax with the frequency 
V (Fig. El^a)) and with the force fmax (Fig. EJ • Inter- 
estingly, Eq. IHl predicts that Axmax varies as g^. This 
prediction is again in excellent agreement with our sim- 
ulation data shown in Fig. El for four different values of 

g- 

Equation |S1 implies that the prefactor C is a mate- 
rial constant that remains independent of all our load- 
ing parameters. Fig. 1201 shows Axmax as a function 
of mg'^ / {fmaxv"^) from different simulations with differ- 
ent values of u, fmax and g. The data are in excellent 
agreement with the linear fit suggested by Eq. |S1 with 
C~0.04. 

The above scaling can be incorporated in the fitting 
form 01 expressing fg as a function of a: — Xj and three 
forces /i, /2 and /s; see Fig. In this fitting form, 
the displacements should be divided by Axmax- We will 
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FIG. 18: Scaling of the maximum displacement /S.Xmax with 
the force amplitude fmax- 



not study here in detail the dependence of /i , /2 and /a 
with respect to loading force parameters mg^ fmin and 
fmax- Our simulations show that /s is independent of 
fmax, but it depends linearly on mg. Theoretically, this 
state corresponds to the limit active state where the ratio 
of principal stresses is a function of the internal angle of 
friction [2j|. On the other hand, the force /2 simply 
scales as fmax and /i depends both on fmin and mg. In 



FIG. 20: Scaling of the maximum displacement /S.Xmax with 
loading parameters from simulations with different values of 
the frequency u (squares), the force amplitude fmax (circles) 
and the gravity g (diamonds). 



our simulations, where /„ 
zero. 



0, the force /i is close to 



VII. COMPACTION RATES 

Equation[3 suggests that the compaction rate rj should 
vary linearly with the frequency v if the total compaction 
per period Api is independent of v. Fig. 1211 shows rj as 
a function of v. We see that only at low frequencies 77 
increases linearly with v. At larger frequencies, beyond 
a characteristic frequency Vc, rj declines with v. The 
largest compaction rate rjmax occurs for u = Vc- This 
implies that, according to Eq. |7| Api is indeed indepen- 
dent of V for V < Vc- The characteristic time Tc = v~^ 
can be interpreted as the minimum time laps required 
for the relaxation of the packing. In fact, in the active 
state, the packing needs a finite rearrangement time Tc 
to achieve a higher level of solid fraction. As long as the 



period t 



is longer than the relaxation time Tc, the 



packing has enough time to relax fully to a more compact 
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FIG. 21: The compaction rate 77 as a function of the frequency 
(circles) fitted by two different functions; see text, v 



state. Then, the compaction Api has its maximum value 
^Pmax- But, if the period r is below Tc, the relaxation 
will be incomplete so that Api < /S.pmax- 

Since the volume change AV^ is proportional to Ax, 
Api follows the same scaling with the frequency as 
the displacement of the retaining wall, i.e. Api oc 



Ap„ 



Hence, from Eq. [7|and imposing the conti- 



nuity at i^ = i^c, we get 



Ap„ 



v 

,,2 






(9) 



This form (labeled 1) is plotted in Fig. l21ltogether with 
the data points. It is remarkable that, although Vc is 
the only fitting parameter, the compaction rate 77 is well 
adjusted by Eq. |21 The prefactor I^pmax/ Pa is — 1-5 x 
10~^, corresponding to ^Pmax — 1-1 x 10^^. 

The arguments behind the proposed form |51 imply a 
sharp transition at v = Vc- This is rather plausible in 
view of the numerical data shown in Fig. 1211 Never- 
theless, it is convenient to construct a single expression 
containing the correct behavior both at low and high fre- 
quencies. The following fitting form provides a good ap- 
proximation as shown also in Fig. [2] (fitting form 2): 



Ap„ 



.-fe-ir 



PO 



1 



(*)■ 



(10) 



We have I'c ~ 6 Hz corresponding to a characteristic 
time Tc = 0.17 s. This time interval is long compared to 
single grain dynamics under gravity. For instance, the 
time required for a small-size grain in our samples to fall 
down a distance equal to its diameter is about 0.002 s. 
Several observations show that collective rearrangements 
in granular media are often a slow process (;27j). Let 
us finally recall that our findings concern the short-time 
behavior {At < 1 s). At longer times, rj declines with 
time, but the scaling with frequency according to Eq. |51 
is expected to hold at each instant of evolution of the 
packing. 



VIII. CONCLUSION 

In this paper, we analyzed the short-time behavior of 
a constrained granular system subjected to vibrational 
dynamics. The vibrations are induced by harmonic vari- 
ation of the force exerted on a free retaining wall be- 
tween zero and a maximum force. The system as a whole 
has a single degree of freedom represented by the hori- 
zontal position of the free wall. This system involves 
a jammed state separating passive (loading) and active 
(unloading) states. The contact dynamics simulations 
were conducted with a rather small number of polygonal 
grains allowing for a systematic study of the dynamics 
and compaction of the material by varying the frequency 
and averaging over several configurations each time. By 
construction, our system is devoid of elastic elements and, 
hence, the behavior is fully governed by collective grain 
rearrangements. 

In the loading phase, the reaction force (exerted by 
the grains on the free wall) rises almost linearly with 
the displacement of the free wall, but it increases con- 
siderably at the end of this phase in transition to the 
jammed state. This force enhancement features the jam- 
ming transition compared to the rest of the passive state. 
The reaction force decreases then in the jammed state, 
balancing thus exactly the driving force, until the lat- 
ter is low enough for the grains to push the free wall 
away under the action of their own weights. This un- 
jamming process occurs smoothly and the reaction force 
increases only slightly but exponentially during the un- 
loading phase. We showed that a rough expression of the 
reaction force as a function of the displacement of the 
free wall with respect to the jamming position, provides 
a good prediction of the dynamics except at the jamming 
and unjamming transients. 

We used dimensional analysis to scale the dynamics 
with the frequency v of oscillations. It was shown that 
the data for frequencies ranging from 1 Hz to 60 Hz, col- 
lapse by scaling the displacements by the inverse square 
of frequency. On the other hand, we studies both numer- 
ically and dimensionally the scaling with loading param- 
eters mg and fmax- 

We also investigated the oscillatory compaction of our 
numerical samples. A small compaction occurs during 
unloading, i.e. during the extension of the sample, fol- 
lowed by a smaller decompaction during loading. The 
compaction rate is nearly constant for short times. It 
was shown that the compaction rate increases linearly 
with frequency up to a characteristic frequency and then 
it declines nearly in inverse proportion to frequency. The 
characteristic frequency was interpreted in terms of the 
time required for the relaxation of a packing in each pe- 
riod to a more compact state by collective grain rear- 
rangements under the action of gravity. The decreasing 
compaction rate as a function of frequency beyond the 
characteristic frequency was explained by arguing that 
only a partial relaxation, inversely proportional to fre- 
quency, could occur at such frequencies. 
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A similar investigation is currently under way with 
polyhedral grains in three dimensions. Our preliminary 
results are consistent with those presented in this paper. 
In view of applications to a wider range of boundary con- 
ditions or diving modes, it is also important to consider 
in detail the characteristic time and the influence of vari- 
ous parameters pertaining to particle properties. Finally, 



long-time behavior and slow evolution of the compaction 
rate may be studied in this framework though more nu- 
merical eff'ort is necessary to reach significant results in 
this case. 

This work was funded by M. Valery from RFF (Roseau 
Ferre Frangais) and the Region Langucdoc-Roussillon 
who are gratefully acknowledged. 
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